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In this paper we consider a system of identical three two-level atoms interacting at resonance 
with a single-mode of the quantized field in a lossless cavity. The initial cavity field is prepared 
in the coherent state while the atoms are taken initially to be either in the uppermost excited 
state "|eee)" or The GHZ-state or the W-state. For this system we investigate different kinds of 
atomic inversion and entanglement, which arise between the different parts of the system due to the 
interaction. Also the relationship, between entanglement and some other nonclassical effects in the 
statistical properties, such as collapses and revivals in the atomic inversion where superharmonic 
effects appear, is discussed. The Q-functions for different cases are discussed. Most remarkably it is 
found that the GHZ-state is more robust against energy losses, showing almost coherent trapping 
and Schrodinger-cat states can not be produced from such state. Also the entanglement of GHZ-state 
is more robust than the W-state. Another interesting feature found is that the state which has no 
pairwise entanglement initially will have a much improvement of such pairwise entanglement through 
the evolution. Sudden death and sudden revival of atoms-pairwise entanglement are produced with 
the W-state. 

PACS numbers: 37.30.-H, 03.67.* 



INTRODUCTION 

The quantum entanglement phenomenon is not only 
one of the most interesting features of the quantum 
theory [ij , that signifies it from the classical theory, but 
also lies at the heart of the new rapidly developing area 
known as the quantum information processing 0- It is 
one of the crucial resources required in the applications 
in this new area of science, which include, quantum 
computation^, quantum teleportation quantum 
dense coding Q and quantum cryptography Q. In 
quantum optics domain, the interaction of an atom with 
a quantized electromagnetic field mode described by the 
Jayens-Cumming model 0, Q leads to an entanglement of 
these two systems such that the total state vector cannot 
be written as a product of the time-dependent atomic 
and field component vectors [§, [13, \lM- To quantify 
entangled states, one should know whether they are pure 
or mixed states. Thus, if the entangled state is in a pure 
state, then it is sufficient to use von Neumann entropy 
as a measure of entanglement. Many efforts have been 
devoted to quantify entanglement, particularly for mixed 
states of a bipartite system, and a number of measures 
have been proposed, such as entanglement of formation, 
relative entropy of entanglement and negativity. The 
Peres-Horodecki criterion for separability [3| leads 
to a natural computable measure of entanglement, 
called negativity [ij, [l^, [l^. It has been proved that 
the negativity Af{p) is an entanglement monotone and 
therefore can be used as a good measure of entanglement 

m. 



On the other hand the cooperative nature of the 
interaction of a quantized radiation field with a system 
of two-level atoms has first been treated by Dicke 
A particular case of the Dicke model, when the atoms 
interact with a single-mode radiation field inside a c avity 
was investigated by Tavis and Cummings TCM [l^, IT9| . 
The interaction between a field and an ensemble of 
atoms develop correlations between the field and the 
atomic systems and between the atomic systems parties 
as well in the course of their dynamics. Quantifying 
this quantum correlations(entanglement) is one of the 
open problems. Multi-qubit systems are of interest to be 
investigated both theoretically and experimentally. In 
recent years, great achievements have been made in the 
application of three-qubit states, because understanding 
of entanglement and dynamics of three-qubit system 
is an important matter, for example, it has many 
applications in quantum cryptography 12011 . quantum 
computation 21 1 and quantum gates 22, 2M- 



In this paper, we consider a system of identical 
three two-level atoms interacting at resonance with a 
single-mode of the quantized field in a lossless cavity. 
The initial cavity field is prepared in the coherent state 
while the atoms may assume different initial states. For 
this system, we investigate some kinds of entanglement 
which arise between the different parts of the system 
due to the interaction. We use the concurrence, the 
generalized I-concurrence and the negativity to quantify 
entanglement between different parties as well as the 
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3-particle residual entanglement. Also we look at the 
relationship between this entanglement and some other 
nonclassical effects in the statistical properties such as 
collapses and revivals in the atomic inversion. The field 
dynamics and atoms-field entanglement are discussed 
through Q-function. 

This paper is organized as follows. In section 2 we 
introduce the system and its solution. Section 3 is de- 
voted to the atomic inversions. In section 4 we study 
the field atom entanglement, two- atom entangle and the 
atomic system purity. Section 5 is devoted to the resid- 
ual 3-particle entanglement. The field dynamics in phase 
space, the possibility of having a cat state and the atoms- 
field entanglement are addressed in section 6. Finally in 
section 7. we conclude the paper with a discussion of the 
results. 



THE MODEL AND ITS TIME EVOLUTION 

We consider a system of three identical two-level atoms 
interacting with a quantized single-mode electromagnetic 
field. Under the rotating-wave approximation and reso- 
nant condition, the Hamiltonian of this system reads: 



2— a,fe,c 



(0 



(1) 



i—a,b,c 



where {h = 1). The terms Hq and Hj represent the 
free and interaction hamiltonians respectively; l^jj is the 
field frequency, and equals the atomic transition fre- 
quency on resonance; a'^_l\ cr'^_^ and cri*^ are the the 
usual raising, lowering and inversion operators for the i*'' 
atom, satisfying [af ^a^l^] = [a^\af] = ±2ai\ 

[(y']l\ cri^^] ~ while a^(a) is the Bose creation (annihila- 
tion) operator for the quantized field mode satisfying the 
commutation relations [a,a^] = 1, and g is the coupling 
constant. Since [Hq^Hj] — 0, it follows that, the Hamil- 
tonian 111) conserves the total number of excitations iV, 
i.e. the total excitation operator 



7V = ata+i 



(2) 



is a constant of motion. This provides a decom- 
position for the system Hilbert space as 7i = 
J2n=o®'^n such that. Ho = {1.9, ff, S; 0)}, Hi = 
{l3,3,ff;l),|5:e,g;0),|5,g,e;0),|e,5r,5;0)} and H2 = 
{|ff,ff,ff;2), \e,g,g;l}, \g,e,g; 1), \g,g,e; 1), |e,e,5;0) 
, |e, (7, e; 0), e, e; 0)} are a one-dimensional , 4- 
dimensional and 7-dimensional cigcnsubspaccs of and 



"^n-i-sl^^o = {\eee;n),\eeg;n+ l),\ege;n+ 1), \gee;n + 
1) , legg; n -h 2) , \geg; n + 2), \gge; n + 2), \ggg; n + 3) } 

are the eight-dimensional eigensubspaces of N. In the 
Hilbert space constituted by the above basis, the interac- 
tion Hi is a diagonal blocks matrix made up of 8 x 8 sub- 
matrices, every sub-matrix represents a subspace corre- 
sponding to a definite excitation number TV. However we 
found that it is more appropriate for this system to use 
the Dicke states as a basis, because of the degener- 
acy produced by the symmetry. For N spin- 1/2 particle 
system, the Dicke states are defined as the states \S,ms) 
that are common eigenstates of both the square of the to- 
tal spin operator and its component along z-axis " the 
quantization axis " with the corresponding eigenval- 
ues iS'(iS' -I- l)h and m^ft. For our system the Dicke states 
in terms of the above mentioned product states are given 

by 

= |eee;n) 

|I?2) = -^ileegin+l) + \ege;n+l) + \gee;n+l)) 

IP3) = -^{\egg;n + 2) + \geg;n + 2) + \gge;n + 2)) 

IPs) = —{2\gee;n+l)~\ege;n+l)-\eeg;n+l)) 

IPs) = --^{\eeg;n+l) - \ege;n+l)) 

\Vj) = -^i\geg;n + 2) + \gge;n + 2) -2\egg;n + 2)) 

\Vs) = -^i\gge;n + 2)-\geg;n + 2)) (3) 

The first four states arc fully symmetric with respect to 
the permutation of particles, while the other four states 
are in fact a two state of mixed symmetry, which corre- 
spond to two degenerate representations, each of dimen- 
sion 2 . By substituting into the Schrodinger equation 



;d|*(r)> 



Hi\^[t)), where r = gt, then the wave func- 
tion generally can be written as 1 4* (t)) = {T)\'Di) . 

We found, as expected, that the equations of the proba- 
bility amplitudes which correspond to the two degener- 
ate states are decoupled from each other and from the 
four amplitudes which correspond to the fully symmetric 
Dicke states. The initial state of the system will be taken 
as follows: 



l*(0)) 

l*a(0)) 



= |vI/„(0))®|vI'/(0)) 

= Ce|e,e,e) -I- -j=Cwi{\e,e,g) + |e,g, e) -f- |.g, e,e)) 



V3 



Cw2{\g,g,e) + \g,e,g) + \e,g,g)) + Cg\g,g,g) 



1^/(0)) = E 



qn\n)] Qn 



,(-|Qo|V2)i^ 



(4) 
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where |Cep + |C^ip + \C^2? + |Cgp = 1. For the above 
initial atomic states the wave function evolve only in- 
side the subspace spanned by the fully symmetric Dicke 
states and the probability amplitudes of the other state 
vanishes, hence our system is equivalent to a four-level 
Dicke atom. This is a generalization to the case of two 
identical atoms [2^, where the antisymmetric states do 
not participate in the dynamics. Solving the coupled four 
equations of the probability amplitudes we get the wave 
function of system at any time 



X(r) 



ri=0 i=l 



UX{0) 



(5) 



where the evolution matrix U is given by 



U2I 1^22 ^23 ^24 
Ud.1 ^32 W33 ^34 
^^41 Ui2 ^43 ^44/ 



(6) 



The explicit form of the matrix elements are given in the 
appendix, where in there we have 



<5 ^ (4/3^+372 + 37/2) 



7 = ^/n+l, 13 = Vn + 2, r] = ^A^+3 (7) 

(n) 

Note that X^ has qn+i-i as a factor, where 1=1,2,3,4. 
By using this wave function, we discuss different aspects 
of the system in what follows. 



ATOMIC INVERSIONS 

Energy is the primary quantity determining the prop- 
erties of physical systems. The atomic inversion and lev- 
els occupation probabilities are the simplest nontrivial 
physical quantities in the atom-field interaction, that dis- 
play the exchange of energy between the field and the 
atoms. More important to our investigation, the atomic 
inversions display the primary nonclassical effect namely, 
the collapse and revival structure \~ 



27, \m\m fror 



which we can have information about the atoms-field en- 
tanglement and disentanglement through the dynamics 
\m, S, [m, [33. Finally the collapses and re- 
vivals will be clearly connected to the evolution of the 
Q-function 3^, 34 1 . In this section we shall discuss these 



quantities. We shall deal with three quantities namely. 
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where, Wxit) is the total atomic inversion, Wa{t) is 
the inversion of the atom a and Pini (t) is the probability 
of occupation of the initial state. The two-atoms inver- 
sion is double the one-atom inversion. In Fig[T]we show 
these quantities for the atoms being initially (a) in the 
excited state |^!n(0)) = |eee), or (b) the atoms initially in 
a GHZ-state 0, |*„(0)) = |GHZ) = -j^Qggg) + |eee)), 
which has the property that tracing over any one 
qubit results in a maximally mixed state containing 
no entanglement between the remaining two qubits or 
(c) in genuine entangled W-state "Werner state" [s^, 
|^'„(0)) = |W) = ^(|ee.g) + \ege) + \gee)), for this state 
when tracing over any one qubit the average remaining 
bipartite entanglement is maximal. Before we dwell in 
the discussion of the figures we present a rough analysis 
about the time dependent quantities in the expressions 
([H [9] and [TO)) . Since we are dealing with a coherent 
state for which n = 100, it is well known that the 
effective excitation will be due to the photons within 
the range |rt — n| < Art = \/E. For these photons and 
those within a reasonably range to them, the radicands 
appearing in the expressions for f3, 7 and 77 in Eq.© 
can be approximated to ^/n and the expressions for /xi, 
^2 are then given by "9n" and "n" respectively. Thus 
the Rabi frequencies appearing in the expressions for 
Uij are to be approximated as 3y/n and y/n respectively. 
Keeping these points in mind, we find the following for 
the case of the atomic initial state |\E'a(0)) = |eee), i.e., 
the three atoms are in their excited states: Regarding 
the quantity Wt{t), the time-dependent summand is 
proportional to j^(cos6Y^r + 15cos2-ynT), containing 
superharmonics which results in three revival times 
^VHj 27rv^ related to the trigonometric 

functions with two equal heights for the 1^* and 2"^^ 
revival and 16 times larger at the 3'''' revival. This 
appears clearly in FiglT](al). On the other hand 
for the single atom population inversion Wa{t), the 
time dependent summand is proportional to cos 2^fnT^ 
which results in a single revival time at 2iry/^ which 
is depicted in Fig[TJ(a2). The PmiXr) has the term 
■^{W + cos6y/nT + 6cos4Y^r + 15cos2y^t) in its 
summand, composed of 2 and 3 superharmoics giving 
rise to revivals times at ^^/E, iry/E, and 2TTy/E. 

The second revival is larger than the I'** and 3'''^ because 
of the coefficients of cos 61/711" and cos 4:^/nT whiles 
the 4*^", at 2tt^/E, is the highest as the four revivals 
coincide. This is shown faithfully in Fig[T](a3) whereas 
the quantity Pmi fluctuates around j^. 
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FIG. 1: Evolution of the total atomic inversions Wt{t), single-atom inversions Wair) and the initial state populations Pi„i{T), 
Eq.(??) for different initial atomic states (a)|eee)-initial state, (b) |GHZ)-initial state and (c) |W)-initial state. The field in a 
coherent state with laP = 100 



Now we look at the initial state |^'a(0)) — 
7f(l.955) + |eee)) Figs[l];ib,2b,3b), and applying 
the same analysis as in the previous case, we find the 
following:- The quantity Wt{t) as well as Wair) is 
almost zero as can be seen from Figlljbl) and Fig[T]Jb2). 
In contrast, the quantity Pini{T) depends on the sum- 
mand of the form |(5 + cos Ay/nr) which results in 
revival times ttv^ and 2tt^/E shown clearly in Fig[TJb3). 
Also it exhibits the fluctuations aroimd |. These results 

o 

show that, starting from the GHZ-state would result in 
coherent trapping and the atoms almost do not interact 
with the field. This is shown in the value of | for the 
probability for the atoms staying in the initial state. 



The case of the atomic initial Werner state namely 



l^'a(O)) 75(1*^^5) + lege) + \gee)) is shown in 

Figs[IJcl,c2,c3). Analysis of the summand of Wt{t) re- 
veals that the time-dependence is of the form (cos Qy/nT+ 
cos2y/riT). This amounts to the revival times 

and 27rv^ with the amplitude at the 3'''' revival 
twice as much as the I''* revival. Whereas the single atom 
population inversion has a single revival at l-n^ffi because 
the time dependent term in the summand is proportional 
to cos2-ynT. It is worthnoting that the fiuctuations in 
the population inversion are the highest for the |eee) ini- 
tial state and the lowest for the GHZ initial state showing 
an almost coherent trapping (37j for the latter in such a 
way to resist exchange of energy with the field, whereas 
the W-state is in the middle, showing a modest degree of 
exchange of the energy with the field. 
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COOPERATIVE AND ATOMS-PAIRWISE 
ENTANGLEMENT 

The system under investigation is a multipartite 
system initially in an over all pure state, occupying the 
following Hilbert space Ti. = Ha ^ Hb He Ti. f with 
dimension 2 (g) 2 (g) 2 (g) oo. Since we consider here only 
the cooperative case, whereas the Hamiltonian is sym- 
metric under atom exchange, thus there are only eight 
nonequivalent partitions of bipartite subsystems namely: 
(i) /(g)(a6c), (/c)® ab, (/c5)(g) a, (ii) f(S){ab), (iii) /(g) a, 
(iv) (fb) g) a, (v) a (g) (be) and (vi) a (g 6, corresponding 
to the field times the whole atomic system, two-atoms 
times the field with the remaining atom, one-atom 
times the field with the other two atoms, the field times 
two-atoms, the field times one-atom, one atom times 
the field with one-atom, one-atom times the other two 
atoms and one-atom times one-atom respectively. 

These partitions except (i), are obtained by tracing 
over one or more Hilbert-subspaces of the total Hilbert 
space and hence they are generally mixed states. For 
convenience and simplicity we will only study the 
entanglement evolution of the partitions (i) and the 
entanglement of the last partition, i.e. atom-atom 
entanglement. In the following section the entanglement 
of the bipartite partition (v), will be discussed together 
with the phenomenon of entanglement sharing 0, [Sq and 
the residual 3-particle entanglement using the negativity. 

For the bipartite partitions in (i), each one of them 
start from an over all pure state and a number of widely 
accepted measures of entanglement are available. An eas- 
ier way to quantify the entanglement in this bipartite par- 
titions is to use the square of the pure-state I-concurrence 
introduced in Ref. 13911 . 



This 
tion 



generalizes 
introduced 



the 

by 



Tr{/5<,ftJ] 
(11) 

original concurrence's no- 
Hill and Wootters ^ 



for pairs of qubits in a joint pure state j^*), to be applied 
for pairs of quantum systems of arbitrary dimension 
di g) ^2 in a joint pure state. The concurrence is defined 
with the help of a superoperator S2, whose action on 
a qubit density operator p = is to flip the spin 

of the qubit density operator S2{p) = p*; p* — |^')(^|, 
where \^) = ay ® cry\^*) , the asterisk denotes the 
complex conjugate and ay is the Pauli matrix. Rungta 



et al. [39| use the formalism for superoperators [41 1 
to generalize the spin-flip superoperator S2 for a qubit 
to a superoperator 5"^ that acts on qudit states (d- 
dimensional states). For defining an I-concurrence, one 
should choose the scaling factor i/^ to be independent 
of d, otherwise, the pure state I-concurrence could be 



changed simply by adding extra, unused dimensions to 
one of the subsystems. To be consistent with the qubit 
concurrence, one should choose = 1. With this choice 
the pure-state I-concurrence runs from zero for product 

states to Imax = \J '^'^™rn'^\ whcrc m = min(di,d2), for a 
maximally entangled state. Henceforth we will use only 
the term I-concurrence when refereing to it. 

On the other hand, Wootters extended the concurrence 
notation to the case of a two qubits in an arbitrary joint 
mixed state, he showed that the entanglement of forma- 
3 3 of an arbitrary two-qubits mixed state p can 



tion 



be written in terms of the minimum average pure-state 
concurrence of ensemble decompositions of p, and he de- 
rived an explicit expression for this minimum in terms of 
the eigenvalues of ppQ. 



C{p) 



c{0, Al - A2 - A3 - A4} 



(12) 



where A^'s are the eigenvalues, in decreasing order, of the 
Hermitian matrix R = ^f^fpp^Jp. Alternatively, one can 
say that the Ai 's are the square roots of the eigenvalues 
of the matrix pp and each Ai is a non-negative real 
number. The spin-flipped state p is obtained by spin 
flipping, namely p = {ay g) ay)p*{ay g) ay), where 
again the asterisk denote the complex conjugate. For a 
pure state j^*), R has only one eigenvalue that may be 
nonzero, namely, C(^') Eq. pTj) . Wootters called this 
minimum average, the concurrence of the mixed state. 
The entanglement of the last partition, i.e. atom-atom 
entanglement can be investigated using this formula. 

It is worthnoting that the remanning partitions are 
all in a mixed state with dimensions 2 g) 00, 4 g) 00 and 
2 (g 4, which can not be quantified by any of the above 
mentioned entanglement measures. Following Tessier et 
al. , we can quantify the entanglement for all the above 
partitions using Osborne's formula [i^ for the I-tangle 
t{pab) for mixed states pab of a pair AB of qudits, of 
dimensions c?^ and ds, having no more than two nonzero 
eigenvalues, i.e., pab with rank no greater than 2. The 
I-tangle t between A and B is given by the expression 

t{pab) = Tt{pabPab} -I- 2 A,„i„(l - TT{p\g}) (13) 

where Xmin is is the smallest eigenvalue of the real 
symmetric 3x3 matrix M as defined in [45| 

In Figl2]we show the I-concurrence of the reduced den- 
sity matrices pabc, Pab and pa, (a) for the |eee) case, (b) 
the I GHZ) case and (c) the |W) case. This type of en- 
tanglement may be called a cooperative entanglement. 
Figs.2(al,a2,a3,bl,cl) begin from zero which corresponds 
to the initial product state while Figs.(b2,b3,c2,c3) start 
from 3-particle maximum entanglement state. They sud- 
denly increase in the collapse region to reach the max- 
imum values, Imax = 1.34,1.23 and 1 for / g) (abc), 
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FIG. 2; The I-concurrence for the three bipartite partitions in (i), (a)]eee)-initial state, (b) |GHZ)-initiaI state and (c) |W)-initiaI 
state. 



(/c)® ab and (fcb)^ a respectively. The exception is the 
entanglement between the field and the total atomic sys- 
tem in the case of the GHZ-statc in FiglSJbl), which did 
not reach this maximum. This reflects the robustness of 
the initial entanglement of this state, i.e., the initial en- 
tanglement between the atoms resists the interaction and 
the energy exchange between them and the field as shown 
in the discussion of the population inversion of last sec- 
tion. When this entanglement becomes rather weak the 
field-atoms entanglement increase rapidly to the maxi- 
mum value Irnax- On the other hand the W-state case 
did not show such robustness. There are many local min- 
ima and maxima asymptotically increasing to reach the 
maximum values, Imax- It is noticed that the GHZ-state 



shows more periodicity than the other two cases. Fi- 
nally, it is clear that the local maxima and minima of the 
f(E) (abc) and (a) (E) fbc figures are opposite to each other. 
FiglSl shows the concurrence of the reduced density ma- 
trix C{pab) for the different initial states. In contrast to 
the above cooperative entanglement, the atom-atom en- 
tanglement asymptotically decreases and vanishes. The 
initially excited atoms and the GHZ case start from zero 
while the |W) case starts from maximum entanglement 



C(pab) = I and then there is a sudden death 4^. 



47, 4 



and sudden rebirth "anabiosis" [49|. This phenomenon 
appears in this case as a bipartite subsystem where the 
third atom and the field are traced. 
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FIG. 3: Evolution of the atoms-pairwise entanglement, by the 
concurrence of Pab For (a)|eee)-initial state, (b) |GHZ)-initial 
state and (c) |W)-initial state. 



THREE-PARTICLE RESIDUAL 
ENTANGLEMENT 

By using the notion of negativity [3j EB], we shall 
study the phenomena of entanglement sharing and the 
residual 3-particle entanglement "three-way entangle- 
ment" 0, [la]. Coffman et al 38 1 introduced the tangle 
notation to quantify the entanglement of three qubits A, 
B, and C in a joint pure state and they discussed en- 
tanglement sharing between three particles. They found 
that unlike classical correlation, quantum entanglement 
cannot be freely shared among the particles. More pre- 
cisely they found that the following inequality 



TA{BC) > TAB + TAC 



(14) 



holds for any three qubits in a joint pure state. Here the 
tangle ta{bc) is the I-concurrence of Eq[Tll between A 
and the other two qubits as one entity, while the other 
two tangles tab and tac are the squares of concurrence 
EqdH between A and B, and A and C respectively. Con- 
sequently they defined the 3-tangle (or residual tangle) 



taibc) = Tab +TAC + tabc 



(15) 



To be read, the entanglement of A with the rest of 
the system is equal to the entanglement of A with B 
alone, plus the entanglement of A with C alone plus the 
3-tangle of the whole system. The 3-tangle quantify the 
residual three-particle entanglement, which cannot be 
accounted for by the pairwise entanglement. For the 
GHZ-state we have 1 =0-1-0-1- tabc- Also they showed 
that the inequality is true in case of three particle 
in a joint mixed state 

On the other hand the Peres-Horodecki criterion for 
separability 13, l3| leads to a natural com put able mea- 
sure of entanglement, called negativity [13, 113,111. The 
negativity is based on the trace norm of the partial trans- 
pose p^-* of the bipartite mixed state pab , and measures 
the degree to which p^^ fails to be positive. The den- 
sity matrices which represent physical systems are non- 
negative matrices with unit trace, the partial transpose 
also satisfies Trlp^^^} = 1, but since it may have negative 
eigenvalues pi" for entangled states" therefore its trace 
norm reads in general 



||i=l + 2|^/i,| = l+A^(p) 



(16) 



where Af{p) is the negativity, i.e. the absolute value of 
twice the sum of the negative eigenvalues. Vidal and 
Werner [13] proved that the negativity A/'(p) is an entan- 
glement monotone and therefore it is a good measure of 
entanglement. Following [l3], the residual three-particle 
entanglement Mabc in our case can be quantified using 
the relation 

^fabc = I^a-bc{Pabc) - I^a-b{Pab) ^ J^a-c{Pac) (17) 



J^a-bciPabc) 
■N'a-b{Pab) 
■N'a^ciPac) 



Pabc 

Pab 
T„ 
Pac 



-1 



-1=1 



Pab 

Tr 

Pac 



-1 



(18) 
(19) 
(20) 



The term Na-bc{Pabc) quantifies the strength of quan- 
tum correlations between the atom "a" and the other 
two atoms. The remanning two terms quantify the pair- 
wise entanglement between the atom " a" and b or c. Note 
that the partial trace operation belongs to the set of local 
operations and classical communication (LOCC) under 
which the entanglement can not increase. Therefore the 
left hand side of Eq.l7 is a residual three-particle entan- 
glement. In our case Afa-b{Pab) = J^a-c{pac) because of 
the symmetry. In FiglH we show the results for the dif- 
ferent initial states. We see that the residual 3-particle 
entanglement changes between local maxima and min- 
ima during the evolution. The local maxima occur at 
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FIG. 4: Evolution of the negativity Ma-bc{pabc), Ma-b{pab) and the residual 3-particles negativity Mabc for different initial 
atomic states (a) |eee)-initial state, (b) |GHZ)-initial state and (c) |W)-initial state. 



half-revival times, when the the atoms became most dis- 
entangled from the field. Asymptotically this residual 
entanglement vanishes and the atomic system becomes 
maximally entangled w\t\i the field. From Figs. [3] and 
|4] it is shown that the negativity measure of two-atom 
entanglement is less than or at most equal to the concur- 
rence as it is conjectured in [l5| . 



Q-FUNCTION 

In this section we return to the entanglement between 
the field and the atomic system but through the field dy- 
namic, namely we will discuss the field quasi- probability 
distribution function, the Q-function [50|, defined as 
Q{p,t) = i(/3|p/(i)|/3). By using Eq.® wc get the 
following form of the Q-function in terms of probability 
amplitudes. 



Q(a,t) 



-|/3| 



2 4 

E 



E 



* Vm+i — 1 



V^nTi^l 



--X. 



in) 



(21) 



Figs l5l61 [7] show the Q-function for different initial 
atomic states at different characteristic times io = 0, 

tl = f\/^, t2 = ^V^, ^3 = TTv/l, = ^V", 

^5 = 2 TTV^. The field is initially in a coherent state with 
average photon number \ao\'^ = n = 100. At r = the 
shape of the Q-function is a Gaussian centered at (-y/H, 0) 
as expected for the initial coherent state. But as time 
develops we found different types of behavior associated 
with each initial atomic state. For the state |eee) 
depicted in FigO we note that the single peak is spfitted 
into two pairs, a pair with large amplitude moving slowly 
in opposite directions and the other pair of the smaller 
amplitudes moving fastly in opposite directions. To shed 
some light on this behavior, rough analysis of formula 
(Plj) . similar to that used in the atomic inversion, shows 
that time dependent factors in the summand behave 
like {cos3(-/n — y/m)T + 3cos(v^ ^ y/Tri)T}. This form 
exhibits clearly the appearance of the pairs and the 
first term represents the faster pair while the second 
term represents the slower pair with amplitude 3 times 
larger than the faster ones. Note the appearance of the 
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to 




FIG. 5: The Q-function for the 

Schrodinger cat-state [s^l, at half-revival (r = t^). The 
case of the state -^{\ggg) + |eee)) is quite interesting. 
Here the peak splits into only two peaks; a small faster 
peak and a large slower peak. The analysis shows that 
the time-dependent factor in the summand behaves as 

|g3i(v77r-VW)r ^ 3gi(v^-v^)r| .^^j^j^,;^ exhibits the twO 

single peaks, they move in opposite directions and the 



tl 




initial atomic state 1*1/(0)) — |eee) 

larger is the slower peak. The third case of the Werner 
state as depicted in Fig[7] shows an opposite behaviour 
to that of the |eee) state. The pair with the larger ampli- 
tude moves faster than the pair with smaller amplitude. 
The time dependent-factor in the summand in this case 
behaves like {3cos3(-y/n — ^/m)T + cos(y^ — y^)t} 
which demonstrates the fast movement of the larger 
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to 




FIG. 6: The Q-function for the initial 

components and the slow movement of the smaller 
components. It is to be remarked that at half revival 
time starting from either the uppermost excited state 
(|eee)) or the Werener state the Schrodinger cat-state 
is produced. On the other hand starting from the 
GHZ-state we find that at half revival time the field 
state returns to its original state. Because as the larger 



tl 




atomic state|»I/(0)) = ■^{\eee) + \ggg)) 

component moves ^ in the phase space the smaller and 
faster component would move ^ on the other direction 
and meets the larger one. The previous behavior is 
also demonstrated in FiglT](b) regarding the atomic 
inversions 

Now looking at the various entanglement evolution 
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to 



tl 





t2 



t3 







FIG. 7: The Q-function for the initial atomic state|'I'(0)) 



%/3 



{\eeg) + \ege) + \gee)) 



Figsl2]and[3]we find that there is a clear connection be- 
tween the Q-function dynamic and various entanglement 
evolutions. For the |eee) and |W) initial states, at half 
the first-revival time ti which correspond to a local min- 



ima of I-Concurrence, / 



/(ahc 



in the leee) case and to the 



minimum value in the |W) case, the faster peaks became 
most apart from each other. Also at that time the resid- 



ual 3-particle entanglement has a local maximum. At the 
first-revival time these two peaks collide and recom- 
bine, and this corresponds to a local minima of If{abc) in 
the two cases, a local minima of the residual 3-particle in 
|eee) case and a maximum value in the |W) case. Also 
we note that the amplitude of these two peaks increase 
significantly in the |W) case. This behavior can be con- 
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nected to Figsl3l(cl),(c2),(c3), where it is clear that at 
that time the entanglement of all the atomic ensem- 
bles with the field is minimum and this is not the case in 
the |eee) state. The last interesting feature we want to 
mention here for the |eee) and Werener cases is that, at 
t5, the Q-function, shows phase squeezing. 

DISCUSSION AND CONCLUSION 

In this paper we consider a system of identical three 
two-level atoms interacting at resonance with a single- 
mode of the quantized field in a lossless cavity. The ini- 
tial cavity field is prepared in the coherent state while 
the atoms are taken in different initial states, namely 
the atoms taken to be in the excited state, |eee), the 
GHZ entangled state and the Werener entangled state. 
For this system wc investigated different kinds of entan- 
glement, atoms-field cooperative and atoms-pairwise en- 
tanglements. We use the concurrence, the generalized 
I-concurrence and the negativity as measures of these 
types of entanglements. The relationship between this 
entanglement and the collapse and revival in the atomic 
inversion is investigated. Also the Q-functions for dif- 
ferent cases are discussed and connected to different en- 
tanglements evolutions of the system. Most noteworthy 
we found that the GHZ-state is more robust against en- 
ergy losses, and showing almost coherent trapping. Also 
one can say that the entanglement of GHZ-state is more 
robust than the W-state. These two different behav- 



iors have been distinctly shown through the study, per- 
fectly depicted in FiglT](bl) and Fig[4l(bl) respectively. 
This suggests that the GHZ-state may show more re- 
sistance to the decoherence phenomena than any other 
three-partite entanglement. Consequently it may be a 
primary candidate for many quantum information tasks, 
which need a three-partite entanglement state. In fact 
the GHZ-state is a resource for many applications, these 
include, quantum secret sharing [51j, open destination 
teleportation [s^ and quantum computation [s^. On 
the other hand for the other two cases (|eee) and W 
states) the Schroinger cat-state is produced. Another in- 
teresting feature is the clear link between the Q-function 
dynamics and various entanglement evolution. Finally 
we found that, while the W-state, is the state with the 
maximal possible bipartite entanglement in the reduced 
two-qubit states, its initial entanglement vanishes very 
rabidly. Moreover the production of such pairwise entan- 
glement through the evolution is very small. This in con- 
trast to the other two cases which have no pairwise entan- 
glement initially but such entanglement increase greatly 
through the evolution. Sudden death and sudden revival 
of atoms-pairwise entanglement are produced with the 
W-state 

APPENDIX : EVOLUTION MATRIX FOR THE 
PROBABILITIES AMPLITUDES 

In this appendix we give the explicit form of the ele- 
ments of the evolution matrix U appearing in Eq. ([6]) 
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^33 
U44 



Ml - 4/32 - 3,72 - 4/32 ^ 3^72 , ^ , 



(Mi - 



(M2 - 



= — i\/3 7 
= 2V3/37 



- M2) 

1 



sin(Y7ir''') + 



M2 ^ 377^ 



Ai2(M2 - Ml) 



sin(Y^/i2 t) 



(Mi - M2) 
1 



(A«2 - Ml) 



cos(yM2 t) 



/T^Kmi - M2) 



sin(Y^r) + 



7^2(^2 - Ml) 



siii(-yM2 t) 



Ml - 377^ s U2~ W I . ^ 

cos(^/^t) + COs(^//ijT) 



(Mi - M2) 



(M2 - Ml) 



/Ml 



(Mi - M2) 
1 



coa{y/jr[T) 



(M2 - Ml) 
1 



(Mi - M2) 
Ml - 372 



cos{y/jr[T) + 

Ml - 37^ 
/m(Mi - M2) 
Ml - 4/32 - 372 



(M2 - Ml) 
M2 - 37' 



sin(VM2T 



cos(^//IJ r) 



(M1-M2) (m2-Mi) 

^ ^ N , M2-372 ■ 



(Mi - M2) 



cos(^/mTt) 



/M2(M2 - Ml) 

M2 - 4/32 - 372 



(M2 - Ml) 



cos(^7IJt) 
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In this paper we consider a system of identical three two-level atoms interacting at resonance 
with a single-mode of the quantized field in a lossless cavity. The initial cavity field is prepared 
in the coherent state while the atoms are taken initially to be either in the uppermost excited 
state "|eee)" or The GHZ-state or the W-state. For this system we investigate different kinds of 
atomic inversion and entanglement, which arise between the different parts of the system due to the 
interaction. Also the relationship, between entanglement and some other nonclassical effects in the 
statistical properties, such as collapses and revivals in the atomic inversion where superharmonic 
effects appear, is discussed. The Q-functions for different cases are discussed. Most remarkably it is 
found that the GHZ-state is more robust against energy losses, showing almost coherent trapping 
and Schrodinger-cat states can not be produced from such state. Also the entanglement of GHZ-state 
is more robust than the W-state. Another interesting feature found is that the state which has no 
pairwise entanglement initially will have a much improvement of such pairwise entanglement through 
the evolution. Sudden death and sudden revival of atoms-pairwise entanglement are produced with 
the W-state. 

PACS numbers: 37.30.-H, 03.67.* 



INTRODUCTION 

The quantum entanglement phenomenon is not only 
one of the most interesting features of the quantum 
theory [ij , that signifies it from the classical theory, but 
also lies at the heart of the new rapidly developing area 
known as the quantum information processing 0- It is 
one of the crucial resources required in the applications 
in this new area of science, which include, quantum 
computation^, quantum teleportation quantum 
dense coding Q and quantum cryptography In 
quantum optics domain, the interaction of an atom with 
a quantized electromagnetic field mode described by the 
Jayens-Cumming model 0, Q leads to an entanglement of 
these two systems such that the total state vector cannot 
be written as a product of the time-dependent atomic 
and field component vectors [§, [13, \lM- To quantify 
entangled states, one should know whether they are pure 
or mixed states. Thus, if the entangled state is in a pure 
state, then it is sufficient to use von Neumann entropy 
as a measure of entanglement. Many efforts have been 
devoted to quantify entanglement, particularly for mixed 
states of a bipartite system, and a number of measures 
have been proposed, such as entanglement of formation, 
relative entropy of entanglement and negativity. The 
Peres-Horodecki criterion for separability [3| leads 
to a natural computable measure of entanglement, 
called negativity [ij, [l^, [l^. It has been proved that 
the negativity Af{p) is an entanglement monotone and 
therefore can be used as a good measure of entanglement 

m. 



On the other hand the cooperative nature of the 
interaction of a quantized radiation field with a system 
of two-level atoms has first been treated by Dicke 
A particular case of the Dicke model, when the atoms 
interact with a single-mode radiation field inside a c avity 
was investigated by Tavis and Cummings TCM [l^, IT9| . 
The interaction between a field and an ensemble of 
atoms develop correlations between the field and the 
atomic systems and between the atomic systems parties 
as well in the course of their dynamics. Quantifying 
this quantum correlations(entanglement) is one of the 
open problems. Multi-qubit systems are of interest to be 
investigated both theoretically and experimentally. In 
recent years, great achievements have been made in the 
application of three-qubit states, because understanding 
of entanglement and dynamics of three-qubit system 
is an important matter, for example, it has many 
applications in quantum cryptography 12011 . quantum 
computation 21 1 and quantum gates 22, 2M- 



In this paper, we consider a system of identical 
three two-level atoms interacting at resonance with a 
single-mode of the quantized field in a lossless cavity. 
The initial cavity field is prepared in the coherent state 
while the atoms may assume different initial states. For 
this system, we investigate some kinds of entanglement 
which arise between the different parts of the system 
due to the interaction. We use the concurrence, the 
generalized I-concurrence and the negativity to quantify 
entanglement between different parties as well as the 
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3-particle residual entanglement. Also we look at the 
relationship between this entanglement and some other 
nonclassical effects in the statistical properties such as 
collapses and revivals in the atomic inversion. The field 
dynamics and atoms-field entanglement are discussed 
through Q-function. 

This paper is organized as follows. In section 2 we 
introduce the system and its solution. Section 3 is de- 
voted to the atomic inversions. In section 4 we study 
the field atom entanglement, two- atom entangle and the 
atomic system purity. Section 5 is devoted to the resid- 
ual 3-particle entanglement. The field dynamics in phase 
space, the possibility of having a cat state and the atoms- 
field entanglement are addressed in section 6. Finally in 
section 7. we conclude the paper with a discussion of the 
results. 



THE MODEL AND ITS TIME EVOLUTION 

We consider a system of three identical two-level atoms 
interacting with a quantized single-mode electromagnetic 
field. Under the rotating-wave approximation and reso- 
nant condition, the Hamiltonian of this system reads: 



2— a,fe,c 



(0 



(1) 



i—a,b,c 



where {h = 1). The terms Hq and Hj represent the 
free and interaction hamiltonians respectively; l^jj is the 
field frequency, and equals the atomic transition fre- 
quency on resonance; a'^_l\ cr'^_^ and cri*^ are the the 
usual raising, lowering and inversion operators for the i*'' 
atom, satisfying [af ^a^l^] = [a^\af] = ±2ai\ 

[(y']l\ cri^^] ~ while a^(a) is the Bose creation (annihila- 
tion) operator for the quantized field mode satisfying the 
commutation relations [a,a^] = 1, and g is the coupling 
constant. Since [Hq^Hj] — 0, it follows that, the Hamil- 
tonian 111) conserves the total number of excitations iV, 
i.e. the total excitation operator 



7V = ata+i 



(2) 



is a constant of motion. This provides a decom- 
position for the system Hilbert space as 7i = 
J2n=o®'^n such that. Ho = {1.9, ff, S; 0)}, Hi = 
{l3,3,ff;l),|5:e,g;0),|5,g,e;0),|e,5r,5;0)} and H2 = 
{|ff,ff,ff;2), \e,g,g;l}, \g,e,g; 1), \g,g,e; 1), |e,e,5;0) 
, |e, (7, e; 0), e, e; 0)} are a one-dimensional , 4- 
dimensional and 7-dimensional cigcnsubspaccs of and 



"^n-i-sl^^o = {\eee;n),\eeg;n+ l),\ege;n+ 1), \gee;n + 
1) , legg; n -h 2) , \geg; n + 2), \gge; n + 2), \ggg; n + 3) } 

are the eight-dimensional eigensubspaces of N. In the 
Hilbert space constituted by the above basis, the interac- 
tion Hi is a diagonal blocks matrix made up of 8 x 8 sub- 
matrices, every sub-matrix represents a subspace corre- 
sponding to a definite excitation number TV. However we 
found that it is more appropriate for this system to use 
the Dicke states as a basis, because of the degener- 
acy produced by the symmetry. For N spin- 1/2 particle 
system, the Dicke states are defined as the states \S,ms) 
that are common eigenstates of both the square of the to- 
tal spin operator and its component along z-axis " the 
quantization axis " with the corresponding eigenval- 
ues iS'(iS' -I- l)h and m^ft. For our system the Dicke states 
in terms of the above mentioned product states are given 

by 

= |eee;n) 

|I?2) = -^ileegin+l) + \ege;n+l) + \gee;n+l)) 

IP3) = -^{\egg;n + 2) + \geg;n + 2) + \gge;n + 2)) 

IPs) = —{2\gee;n+l)~\ege;n+l)-\eeg;n+l)) 

IPs) = --^{\eeg;n+l) - \ege;n+l)) 

\Vj) = -^i\geg;n + 2) + \gge;n + 2) -2\egg;n + 2)) 

\Vs) = -^i\gge;n + 2)-\geg;n + 2)) (3) 

The first four states arc fully symmetric with respect to 
the permutation of particles, while the other four states 
are in fact a two state of mixed symmetry, which corre- 
spond to two degenerate representations, each of dimen- 
sion 2 . By substituting into the Schrodinger equation 



;d|*(r)> 



Hi\^[t)), where r = gt, then the wave func- 
tion generally can be written as 1 4* (t)) = {T)\'Di) . 

We found, as expected, that the equations of the proba- 
bility amplitudes which correspond to the two degener- 
ate states are decoupled from each other and from the 
four amplitudes which correspond to the fully symmetric 
Dicke states. The initial state of the system will be taken 
as follows: 



l*(0)) 

l*a(0)) 



= |vI/„(0))®|vI'/(0)) 

= Ce|e,e,e) -I- -j=Cwi{\e,e,g) + |e,g, e) -f- |.g, e,e)) 



V3 



Cw2{\g,g,e) + \g,e,g) + \e,g,g)) + Cg\g,g,g) 



1^/(0)) = E 



qn\n)] Qn 



,(-|Qo|V2)i^ 



(4) 
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where |Cep + |C^ip + \C^2? + |Cgp = 1. For the above 
initial atomic states the wave function evolve only in- 
side the subspace spanned by the fully symmetric Dicke 
states and the probability amplitudes of the other state 
vanishes, hence our system is equivalent to a four-level 
Dicke atom. This is a generalization to the case of two 
identical atoms [2^, where the antisymmetric states do 
not participate in the dynamics. Solving the coupled four 
equations of the probability amplitudes we get the wave 
function of system at any time 



X(r) 



ri=0 i=l 



UX{0) 



(5) 



where the evolution matrix U is given by 



U2I 1^22 ^23 ^24 
Ud.1 ^32 W33 ^34 
^^41 Ui2 ^43 ^44/ 



(6) 



The explicit form of the matrix elements are given in the 
appendix, where in there we have 



<5 ^ (4/3^+372 + 37/2) 



7 = ^/n+l, 13 = Vn + 2, r] = ^A^+3 (7) 

(n) 

Note that X^ has qn+i-i as a factor, where 1=1,2,3,4. 
By using this wave function, we discuss different aspects 
of the system in what follows. 



ATOMIC INVERSIONS 

Energy is the primary quantity determining the prop- 
erties of physical systems. The atomic inversion and lev- 
els occupation probabilities are the simplest nontrivial 
physical quantities in the atom-field interaction, that dis- 
play the exchange of energy between the field and the 
atoms. More important to our investigation, the atomic 
inversions display the primary nonclassical effect namely, 
the collapse and revival structure \~ 



27, \m\m fror 



which we can have information about the atoms-field en- 
tanglement and disentanglement through the dynamics 
\m, S, [m, [33. Finally the collapses and re- 
vivals will be clearly connected to the evolution of the 
Q-function 3^, 34 1 . In this section we shall discuss these 



quantities. We shall deal with three quantities namely. 
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IX, 



(n),2_l|^(n)|2} 



(9) 



(10) 



where, Wxit) is the total atomic inversion, Wa{t) is 
the inversion of the atom a and Pini (t) is the probability 
of occupation of the initial state. The two-atoms inver- 
sion is double the one-atom inversion. In Fig[T]we show 
these quantities for the atoms being initially (a) in the 
excited state |^!n(0)) = |eee), or (b) the atoms initially in 
a GHZ-state 0, |*„(0)) = |GHZ) = -j^Qggg) + |eee)), 
which has the property that tracing over any one 
qubit results in a maximally mixed state containing 
no entanglement between the remaining two qubits or 
(c) in genuine entangled W-state "Werner state" [s^, 
|^'„(0)) = |W) = ^(|ee.g) + \ege) + \gee)), for this state 
when tracing over any one qubit the average remaining 
bipartite entanglement is maximal. Before we dwell in 
the discussion of the figures we present a rough analysis 
about the time dependent quantities in the expressions 
([H [9] and [TO)) . Since we are dealing with a coherent 
state for which n = 100, it is well known that the 
effective excitation will be due to the photons within 
the range |rt — n| < Art = \/E. For these photons and 
those within a reasonably range to them, the radicands 
appearing in the expressions for f3, 7 and 77 in Eq.© 
can be approximated to ^/n and the expressions for /xi, 
^2 are then given by "9n" and "n" respectively. Thus 
the Rabi frequencies appearing in the expressions for 
Uij are to be approximated as 3y/n and y/n respectively. 
Keeping these points in mind, we find the following for 
the case of the atomic initial state |\E'a(0)) = |eee), i.e., 
the three atoms are in their excited states: Regarding 
the quantity Wt{t), the time-dependent summand is 
proportional to j^(cos6Y^r + 15cos2-ynT), containing 
superharmonics which results in three revival times 
^VHj 27rv^ related to the trigonometric 

functions with two equal heights for the 1^* and 2"^^ 
revival and 16 times larger at the 3'''' revival. This 
appears clearly in FiglT](al). On the other hand 
for the single atom population inversion Wa{t), the 
time dependent summand is proportional to cos 2^fnT^ 
which results in a single revival time at 2iry/^ which 
is depicted in Fig[TJ(a2). The PmiXr) has the term 
■^{W + cos6y/nT + 6cos4Y^r + 15cos2y^t) in its 
summand, composed of 2 and 3 superharmoics giving 
rise to revivals times at ^^/E, iry/E, and 2TTy/E. 

The second revival is larger than the I'** and 3'''^ because 
of the coefficients of cos 61/711" and cos 4:^/nT whiles 
the 4*^", at 2tt^/E, is the highest as the four revivals 
coincide. This is shown faithfully in Fig[T](a3) whereas 
the quantity Pmi fluctuates around j^. 
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FIG. 1: Evolution of the total atomic inversions Wt{t), single-atom inversions Wair) and the initial state populations Pi„i{T), 
Eq.(??) for different initial atomic states (a)|eee)-initial state, (b) |GHZ)-initial state and (c) |W)-initial state. The field in a 
coherent state with laP = 100 



Now we look at the initial state |^'a(0)) — 
7f(l.955) + |eee)) Figs[l];ib,2b,3b), and applying 
the same analysis as in the previous case, we find the 
following:- The quantity Wt{t) as well as Wair) is 
almost zero as can be seen from Figlljbl) and Fig[T]Jb2). 
In contrast, the quantity Pini{T) depends on the sum- 
mand of the form |(5 + cos Ay/nr) which results in 
revival times ttv^ and 2tt^/E shown clearly in Fig[TJb3). 
Also it exhibits the fluctuations aroimd |. These results 

o 

show that, starting from the GHZ-state would result in 
coherent trapping and the atoms almost do not interact 
with the field. This is shown in the value of | for the 
probability for the atoms staying in the initial state. 



The case of the atomic initial Werner state namely 



l^'a(O)) 75(1*^^5) + lege) + \gee)) is shown in 

Figs[IJcl,c2,c3). Analysis of the summand of Wt{t) re- 
veals that the time-dependence is of the form (cos Qy/nT+ 
cos2y/riT). This amounts to the revival times 

and 27rv^ with the amplitude at the 3'''' revival 
twice as much as the I''* revival. Whereas the single atom 
population inversion has a single revival at l-n^ffi because 
the time dependent term in the summand is proportional 
to cos2-ynT. It is worthnoting that the fiuctuations in 
the population inversion are the highest for the |eee) ini- 
tial state and the lowest for the GHZ initial state showing 
an almost coherent trapping (37j for the latter in such a 
way to resist exchange of energy with the field, whereas 
the W-state is in the middle, showing a modest degree of 
exchange of the energy with the field. 
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COOPERATIVE AND ATOMS-PAIRWISE 
ENTANGLEMENT 

The system under investigation is a multipartite 
system initially in an over all pure state, occupying the 
following Hilbert space Ti. = Ha ^ Hb He Ti. f with 
dimension 2 (g) 2 (g) 2 (g) oo. Since we consider here only 
the cooperative case, whereas the Hamiltonian is sym- 
metric under atom exchange, thus there are only eight 
nonequivalent partitions of bipartite subsystems namely: 
(i) /(g)(a6c), (/c)® ab, (/c5)(g) a, (ii) f(S){ab), (iii) /(g) a, 
(iv) (fb) g) a, (v) a (g) (be) and (vi) a (g 6, corresponding 
to the field times the whole atomic system, two-atoms 
times the field with the remaining atom, one-atom 
times the field with the other two atoms, the field times 
two-atoms, the field times one-atom, one atom times 
the field with one-atom, one-atom times the other two 
atoms and one-atom times one-atom respectively. 

These partitions except (i), are obtained by tracing 
over one or more Hilbert-subspaces of the total Hilbert 
space and hence they are generally mixed states. For 
convenience and simplicity we will only study the 
entanglement evolution of the partitions (i) and the 
entanglement of the last partition, i.e. atom-atom 
entanglement. In the following section the entanglement 
of the bipartite partition (v), will be discussed together 
with the phenomenon of entanglement sharing 0, [Sq and 
the residual 3-particle entanglement using the negativity. 

For the bipartite partitions in (i), each one of them 
start from an over all pure state and a number of widely 
accepted measures of entanglement are available. An eas- 
ier way to quantify the entanglement in this bipartite par- 
titions is to use the square of the pure-state I-concurrence 
introduced in Ref. 13911 . 



This 
tion 



generalizes 
introduced 



the 

by 



Tr{/5<,ftJ] 
(11) 

original concurrence's no- 
Hill and Wootters ^ 



for pairs of qubits in a joint pure state j^*), to be applied 
for pairs of quantum systems of arbitrary dimension 
di g) ^2 in a joint pure state. The concurrence is defined 
with the help of a superoperator S2, whose action on 
a qubit density operator p = is to flip the spin 

of the qubit density operator S2{p) = p*; p* — |^')(^|, 
where \^) = ay ® cry\^*) , the asterisk denotes the 
complex conjugate and ay is the Pauli matrix. Rungta 



et al. [39| use the formalism for superoperators [41 1 
to generalize the spin-flip superoperator S2 for a qubit 
to a superoperator 5"^ that acts on qudit states (d- 
dimensional states). For defining an I-concurrence, one 
should choose the scaling factor i/^ to be independent 
of d, otherwise, the pure state I-concurrence could be 



changed simply by adding extra, unused dimensions to 
one of the subsystems. To be consistent with the qubit 
concurrence, one should choose = 1. With this choice 
the pure-state I-concurrence runs from zero for product 

states to Imax = \J '^'^™rn'^\ whcrc m = min(di,d2), for a 
maximally entangled state. Henceforth we will use only 
the term I-concurrence when refereing to it. 

On the other hand, Wootters extended the concurrence 
notation to the case of a two qubits in an arbitrary joint 
mixed state, he showed that the entanglement of forma- 
3 3 of an arbitrary two-qubits mixed state p can 



tion 



be written in terms of the minimum average pure-state 
concurrence of ensemble decompositions of p, and he de- 
rived an explicit expression for this minimum in terms of 
the eigenvalues of ppQ. 



C{p) 



c{0, Al - A2 - A3 - A4} 



(12) 



where A^'s are the eigenvalues, in decreasing order, of the 
Hermitian matrix R = ^f^fpp^Jp. Alternatively, one can 
say that the Ai 's are the square roots of the eigenvalues 
of the matrix pp and each Ai is a non-negative real 
number. The spin-flipped state p is obtained by spin 
flipping, namely p = {ay g) ay)p*{ay g) ay), where 
again the asterisk denote the complex conjugate. For a 
pure state j^*), R has only one eigenvalue that may be 
nonzero, namely, C(^') Eq. pTj) . Wootters called this 
minimum average, the concurrence of the mixed state. 
The entanglement of the last partition, i.e. atom-atom 
entanglement can be investigated using this formula. 

It is worthnoting that the remanning partitions are 
all in a mixed state with dimensions 2 g) 00, 4 g) 00 and 
2 (g 4, which can not be quantified by any of the above 
mentioned entanglement measures. Following Tessier et 
al. , we can quantify the entanglement for all the above 
partitions using Osborne's formula [i^ for the I-tangle 
t{pab) for mixed states pab of a pair AB of qudits, of 
dimensions c?^ and ds, having no more than two nonzero 
eigenvalues, i.e., pab with rank no greater than 2. The 
I-tangle t between A and B is given by the expression 

t{pab) = Tt{pabPab} -I- 2 A,„i„(l - TT{p\g}) (13) 

where Xmin is is the smallest eigenvalue of the real 
symmetric 3x3 matrix M as defined in [45| 

In Figl2]we show the I-concurrence of the reduced den- 
sity matrices pabc, Pab and pa, (a) for the |eee) case, (b) 
the I GHZ) case and (c) the |W) case. This type of en- 
tanglement may be called a cooperative entanglement. 
Figs.2(al,a2,a3,bl,cl) begin from zero which corresponds 
to the initial product state while Figs.(b2,b3,c2,c3) start 
from 3-particle maximum entanglement state. They sud- 
denly increase in the collapse region to reach the max- 
imum values, Imax = 1.34,1.23 and 1 for / g) (abc), 
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FIG. 2; The I-concurrence for the three bipartite partitions in (i), (a)]eee)-initial state, (b) |GHZ)-initiaI state and (c) |W)-initiaI 
state. 



(/c)® ab and (fcb)^ a respectively. The exception is the 
entanglement between the field and the total atomic sys- 
tem in the case of the GHZ-statc in FiglSJbl), which did 
not reach this maximum. This reflects the robustness of 
the initial entanglement of this state, i.e., the initial en- 
tanglement between the atoms resists the interaction and 
the energy exchange between them and the field as shown 
in the discussion of the population inversion of last sec- 
tion. When this entanglement becomes rather weak the 
field-atoms entanglement increase rapidly to the maxi- 
mum value Irnax- On the other hand the W-state case 
did not show such robustness. There are many local min- 
ima and maxima asymptotically increasing to reach the 
maximum values, Imax- It is noticed that the GHZ-state 



shows more periodicity than the other two cases. Fi- 
nally, it is clear that the local maxima and minima of the 
f(E) (abc) and (a) (E) fbc figures are opposite to each other. 
FiglSl shows the concurrence of the reduced density ma- 
trix C{pab) for the different initial states. In contrast to 
the above cooperative entanglement, the atom-atom en- 
tanglement asymptotically decreases and vanishes. The 
initially excited atoms and the GHZ case start from zero 
while the |W) case starts from maximum entanglement 



C(pab) = I and then there is a sudden death 4^. 



47, 4 



and sudden rebirth "anabiosis" [49|. This phenomenon 
appears in this case as a bipartite subsystem where the 
third atom and the field are traced. 
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FIG. 3: Evolution of the atoms-pairwise entanglement, by the 
concurrence of Pab For (a)|eee)-initial state, (b) |GHZ)-initial 
state and (c) |W)-initial state. 



THREE-PARTICLE RESIDUAL 
ENTANGLEMENT 

By using the notion of negativity [3j EB], we shall 
study the phenomena of entanglement sharing and the 
residual 3-particle entanglement "three-way entangle- 
ment" 0, [la]. Coffman et al 38 1 introduced the tangle 
notation to quantify the entanglement of three qubits A, 
B, and C in a joint pure state and they discussed en- 
tanglement sharing between three particles. They found 
that unlike classical correlation, quantum entanglement 
cannot be freely shared among the particles. More pre- 
cisely they found that the following inequality 



TA{BC) > TAB + TAC 



(14) 



holds for any three qubits in a joint pure state. Here the 
tangle ta{bc) is the I-concurrence of Eq[Tll between A 
and the other two qubits as one entity, while the other 
two tangles tab and tac are the squares of concurrence 
EqdH between A and B, and A and C respectively. Con- 
sequently they defined the 3-tangle (or residual tangle) 



taibc) = Tab +TAC + tabc 



(15) 



To be read, the entanglement of A with the rest of 
the system is equal to the entanglement of A with B 
alone, plus the entanglement of A with C alone plus the 
3-tangle of the whole system. The 3-tangle quantify the 
residual three-particle entanglement, which cannot be 
accounted for by the pairwise entanglement. For the 
GHZ-state we have 1 =0-1-0-1- tabc- Also they showed 
that the inequality is true in case of three particle 
in a joint mixed state 

On the other hand the Peres-Horodecki criterion for 
separability 13, l3| leads to a natural com put able mea- 
sure of entanglement, called negativity [13, 113,111. The 
negativity is based on the trace norm of the partial trans- 
pose p^-* of the bipartite mixed state pab , and measures 
the degree to which p^^ fails to be positive. The den- 
sity matrices which represent physical systems are non- 
negative matrices with unit trace, the partial transpose 
also satisfies Trlp^^^} = 1, but since it may have negative 
eigenvalues pi" for entangled states" therefore its trace 
norm reads in general 



||i=l + 2|^/i,| = l+A^(p) 



(16) 



where Af{p) is the negativity, i.e. the absolute value of 
twice the sum of the negative eigenvalues. Vidal and 
Werner [13] proved that the negativity A/'(p) is an entan- 
glement monotone and therefore it is a good measure of 
entanglement. Following [l3], the residual three-particle 
entanglement Mabc in our case can be quantified using 
the relation 

^fabc = I^a-bc{Pabc) - I^a-b{Pab) ^ J^a-c{Pac) (17) 



J^a-bciPabc) 
■N'a-b{Pab) 
■N'a^ciPac) 



Pabc 

Pab 
T„ 
Pac 



-1 



-1=1 



Pab 

Tr 

Pac 



-1 



(18) 
(19) 
(20) 



The term Na-bc{Pabc) quantifies the strength of quan- 
tum correlations between the atom "a" and the other 
two atoms. The remanning two terms quantify the pair- 
wise entanglement between the atom " a" and b or c. Note 
that the partial trace operation belongs to the set of local 
operations and classical communication (LOCC) under 
which the entanglement can not increase. Therefore the 
left hand side of Eq.l7 is a residual three-particle entan- 
glement. In our case Afa-b{Pab) = J^a-c{pac) because of 
the symmetry. In FiglH we show the results for the dif- 
ferent initial states. We see that the residual 3-particle 
entanglement changes between local maxima and min- 
ima during the evolution. The local maxima occur at 
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FIG. 4: Evolution of the negativity Ma-bc{pabc), Ma-b{pab) and the residual 3-particles negativity Mabc for different initial 
atomic states (a) |eee)-initial state, (b) |GHZ)-initial state and (c) |W)-initial state. 



half-revival times, when the the atoms became most dis- 
entangled from the field. Asymptotically this residual 
entanglement vanishes and the atomic system becomes 
maximally entangled w\t\i the field. From Figs. [3] and 
|4] it is shown that the negativity measure of two-atom 
entanglement is less than or at most equal to the concur- 
rence as it is conjectured in [l5| . 



Q-FUNCTION 

In this section we return to the entanglement between 
the field and the atomic system but through the field dy- 
namic, namely we will discuss the field quasi- probability 
distribution function, the Q-function [50|, defined as 
Q{p,t) = i(/3|p/(i)|/3). By using Eq.® wc get the 
following form of the Q-function in terms of probability 
amplitudes. 



Q(a,t) 



-|/3| 



2 4 

E 



E 



* Vm+i — 1 



V^nTi^l 



--X. 



in) 



(21) 



Figs l5l61 [7] show the Q-function for different initial 
atomic states at different characteristic times io = 0, 

tl = f\/^, t2 = ^V^, ^3 = TTv/l, = ^V", 

^5 = 2 TTV^. The field is initially in a coherent state with 
average photon number \ao\'^ = n = 100. At r = the 
shape of the Q-function is a Gaussian centered at (-y/H, 0) 
as expected for the initial coherent state. But as time 
develops we found different types of behavior associated 
with each initial atomic state. For the state |eee) 
depicted in FigO we note that the single peak is spfitted 
into two pairs, a pair with large amplitude moving slowly 
in opposite directions and the other pair of the smaller 
amplitudes moving fastly in opposite directions. To shed 
some light on this behavior, rough analysis of formula 
(Plj) . similar to that used in the atomic inversion, shows 
that time dependent factors in the summand behave 
like {cos3(-/n — y/m)T + 3cos(v^ ^ y/Tri)T}. This form 
exhibits clearly the appearance of the pairs and the 
first term represents the faster pair while the second 
term represents the slower pair with amplitude 3 times 
larger than the faster ones. Note the appearance of the 
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to 




FIG. 5: The Q-function for the 

Schrodinger cat-state [s^l, at half-revival (r = t^). The 
case of the state -^{\ggg) + |eee)) is quite interesting. 
Here the peak splits into only two peaks; a small faster 
peak and a large slower peak. The analysis shows that 
the time-dependent factor in the summand behaves as 

|g3i(v77r-VW)r ^ 3gi(v^-v^)r| .^^j^j^,;^ exhibits the twO 

single peaks, they move in opposite directions and the 



tl 




initial atomic state 1*1/(0)) — |eee) 

larger is the slower peak. The third case of the Werner 
state as depicted in Fig[7] shows an opposite behaviour 
to that of the |eee) state. The pair with the larger ampli- 
tude moves faster than the pair with smaller amplitude. 
The time dependent-factor in the summand in this case 
behaves like {3cos3(-y/n — ^/m)T + cos(y^ — y^)t} 
which demonstrates the fast movement of the larger 
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to 




FIG. 6: The Q-function for the initial 

components and the slow movement of the smaller 
components. It is to be remarked that at half revival 
time starting from either the uppermost excited state 
(|eee)) or the Werener state the Schrodinger cat-state 
is produced. On the other hand starting from the 
GHZ-state we find that at half revival time the field 
state returns to its original state. Because as the larger 



tl 




atomic state|»I/(0)) = ■^{\eee) + \ggg)) 

component moves ^ in the phase space the smaller and 
faster component would move ^ on the other direction 
and meets the larger one. The previous behavior is 
also demonstrated in FiglT](b) regarding the atomic 
inversions 

Now looking at the various entanglement evolution 
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tl 





t2 



t3 







FIG. 7: The Q-function for the initial atomic state|'I'(0)) 



%/3 



{\eeg) + \ege) + \gee)) 



Figsl2]and[3]we find that there is a clear connection be- 
tween the Q-function dynamic and various entanglement 
evolutions. For the |eee) and |W) initial states, at half 
the first-revival time ti which correspond to a local min- 



ima of I-Concurrence, / 



/(ahc 



in the leee) case and to the 



minimum value in the |W) case, the faster peaks became 
most apart from each other. Also at that time the resid- 



ual 3-particle entanglement has a local maximum. At the 
first-revival time these two peaks collide and recom- 
bine, and this corresponds to a local minima of If{abc) in 
the two cases, a local minima of the residual 3-particle in 
|eee) case and a maximum value in the |W) case. Also 
we note that the amplitude of these two peaks increase 
significantly in the |W) case. This behavior can be con- 
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nected to Figsl3l(cl),(c2),(c3), where it is clear that at 
that time the entanglement of all the atomic ensem- 
bles with the field is minimum and this is not the case in 
the |eee) state. The last interesting feature we want to 
mention here for the |eee) and Werener cases is that, at 
t5, the Q-function, shows phase squeezing. 

DISCUSSION AND CONCLUSION 

In this paper we consider a system of identical three 
two-level atoms interacting at resonance with a single- 
mode of the quantized field in a lossless cavity. The ini- 
tial cavity field is prepared in the coherent state while 
the atoms are taken in different initial states, namely 
the atoms taken to be in the excited state, |eee), the 
GHZ entangled state and the Werener entangled state. 
For this system wc investigated different kinds of entan- 
glement, atoms-field cooperative and atoms-pairwise en- 
tanglements. We use the concurrence, the generalized 
I-concurrence and the negativity as measures of these 
types of entanglements. The relationship between this 
entanglement and the collapse and revival in the atomic 
inversion is investigated. Also the Q-functions for dif- 
ferent cases are discussed and connected to different en- 
tanglements evolutions of the system. Most noteworthy 
we found that the GHZ-state is more robust against en- 
ergy losses, and showing almost coherent trapping. Also 
one can say that the entanglement of GHZ-state is more 
robust than the W-state. These two different behav- 



iors have been distinctly shown through the study, per- 
fectly depicted in FiglT](bl) and Fig[4l(bl) respectively. 
This suggests that the GHZ-state may show more re- 
sistance to the decoherence phenomena than any other 
three-partite entanglement. Consequently it may be a 
primary candidate for many quantum information tasks, 
which need a three-partite entanglement state. In fact 
the GHZ-state is a resource for many applications, these 
include, quantum secret sharing [51j, open destination 
teleportation [s^ and quantum computation [s^. On 
the other hand for the other two cases (|eee) and W 
states) the Schroinger cat-state is produced. Another in- 
teresting feature is the clear link between the Q-function 
dynamics and various entanglement evolution. Finally 
we found that, while the W-state, is the state with the 
maximal possible bipartite entanglement in the reduced 
two-qubit states, its initial entanglement vanishes very 
rabidly. Moreover the production of such pairwise entan- 
glement through the evolution is very small. This in con- 
trast to the other two cases which have no pairwise entan- 
glement initially but such entanglement increase greatly 
through the evolution. Sudden death and sudden revival 
of atoms-pairwise entanglement are produced with the 
W-state 

APPENDIX : EVOLUTION MATRIX FOR THE 
PROBABILITIES AMPLITUDES 

In this appendix we give the explicit form of the ele- 
ments of the evolution matrix U appearing in Eq. ([6]) 
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^33 
U44 



Ml - 4/32 - 3,72 - 4/32 ^ 3^72 , ^ , 



(Mi - 



(M2 - 



= — i\/3 7 
= 2V3/37 



- M2) 

1 



sin(Y7ir''') + 



M2 ^ 377^ 



Ai2(M2 - Ml) 



sin(Y^/i2 t) 



(Mi - M2) 
1 



(A«2 - Ml) 



cos(yM2 t) 



/T^Kmi - M2) 



sin(Y^r) + 



7^2(^2 - Ml) 



siii(-yM2 t) 



Ml - 377^ s U2~ W I . ^ 

cos(^/^t) + COs(^//ijT) 



(Mi - M2) 



(M2 - Ml) 



/Ml 



(Mi - M2) 
1 



coa{y/jr[T) 



(M2 - Ml) 
1 



(Mi - M2) 
Ml - 372 



cos{y/jr[T) + 

Ml - 37^ 
/m(Mi - M2) 
Ml - 4/32 - 372 



(M2 - Ml) 
M2 - 37' 



sin(VM2T 



cos(^//IJ r) 



(M1-M2) (m2-Mi) 

^ ^ N , M2-372 ■ 



(Mi - M2) 



cos(^/mTt) 



/M2(M2 - Ml) 

M2 - 4/32 - 372 



(M2 - Ml) 



cos(^7IJt) 
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